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Abstract 

We study flavor changing processes T — > B/BX S and J/^ — > D/DX U in the B factories and the 
Tau-Charm factories. In the standard model, these processes are predicted to be unobservable, so they 
serve as a probe of the new physics. We first perform a model independent analysis, then examine 
the predictions of models; such as TopColor models, MSSM with R-parity violation and the two Higgs 
doublet model; for the branching ratios of T — * B/BX S and J/ip — > D/DX U . We find that these 
branching ratios could be as large as 10 -6 and 1CP 5 in the presence of new physics. 



1 Introduction 



The possibility of observing large CP violating asymmetries in the decay of B mesons has motivated 
the construction of high luminosity B factories at several of the world's high energy physics laboratories. 
These B factories will be producing roughly about 10 8 Upsilons. Meanwhile BES has already accumulated 
9 x 10 6 J/ip and plans to increase the number to 5 x 10 7 in the near future. An interesting question, that 
we investigate in this paper, is whether the large sample of the T and the J ftp can be used to probe flavor 
changing processes in the decays of T and J ftp- I n particular we look at the flavor changing processes 
T —* B/BX S and J/ip — » D/DX U , from the underlying b —* s and c — > u quark transitions. For the 
quarkonium system, these flavor changing processes are expected to be much smaller than in the case of 
decays of the B or the D meson because of the larger decay widths of the bottonium and the charmonium 
systems which decay via the strong interactions. Indeed the standard model contributions to T — ► B/BX S 
and J ftp — > D/DX U are tiny. However, new physics may enhance the branching ratios for these processes. 
Whether this enhancement maybe sufficient for these processes to be observable in the next round of 
experiments is the subject of this work. Invisible decays of T and J/ip resonances in the standard model 
and beyond have been studied recently [|l|]. 

Non leptonic decays of heavy quarkonium systems can be more reliably calculated than the non leptonic 
decays of the heavy mesons. A consistent and systematic formalism to handle heavy quarkonium decays 
is available in NRQCD Q which is missing for the heavy mesons. As in the meson system Q it is more 
fruitful to concentrate on quasi- inclusive processes like T — * B/BX S and J ftp — > D/DX U because they can 
be calculated with less theoretical uncertainty and have larger branching ratios than the purely exclusive 
quarkonium non leptonic decays. The branching ratios of exclusive flavor changing non leptonic decays of 
T and J ftp in the standard model have been calculated and found to be very small ||]. 

We begin with a model independent description of the processes T — ► B/BX S and J/ip — > D/DX U . 
In the standard model these decays can proceed through tree and penguin processes. For new physics 
contribution to these processes we concentrate on four quark operators of the type ~sbbb and uccc. We 
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choose the currents in the four quark operators to be scalars and so these operators may arise through 
the exchange of a heavy scalar for e.g a Higgs or a leptoquark in some model of new physics. These four 
quark operators, at the one loop level, generate effective s&{g,7, Z} and uc{g,^, Z} vertices which would 
effect the flavor changing decays of the B and the D mesons. The effective vertices for an on shell g and 
7 vanish and so there is no contribution to b — > 57 or c — > wy. We can however put constraints on these 
operators by considering the processes b — ► sl + l~ and c — ► ul + l~ . The constrained operators can then be 
used to calculate the branching ratios for T — > B/BX S and J /if} — > D/DX U . 

We then consider some models that may generate the kind of four quark operators described above. 
A few examples of models where these operators can be generated are top color models, MSSM with R 
parity violation and a general two Higgs doublet model without any discrete symmetry. In some cases 
constraints on the parameters that appear in the prediction for the branching ratios for T — ► B/BX S and 
J/tjj — > D/DX U are already available. In other cases the parameters are constrained, as in our model 
independent analysis, from the processes b — > sl + l~ and c — ► ul + l~ . 

In the sections which follow, we describe the effective Hamiltonian for the T — ► B/BX S and J/ip — > 
D/DX U . Next we describe the calculation of the matrix elements and decay rates for these processes. We 
then discuss the calculation of the effective sb{g,^, Z} and uc{g,^,Z} vertices and constraints from the 
processes b —* sl + l~ and c — > ul + l~ . This is followed by a description of some models that can generate 
the new four quark operators in the effective Hamiltonian for T — > B/BX S and J ftp — » D/DX U . Finally 
we present our results and conclusions. 

2 Effective Hamiltonian 

In this section we present the effective Hamiltonian for T decays. The effective Hamiltonian for charmonium 
decays can be written down by making obvious changes. In the Standard Model (SM) the amplitudes for 

hadronic T decays of the type bb — > sb + sb are generated by the following effective Hamiltonian ||, |6| : 

f-i 10 
H eff = ^iVfbV^Olf+^O^-^VubV^ + V^V^ (1) 

v 1 i=3 
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where the superscript u, c, t indicates the internal quark, / can be u or c quark, q can be either a d or a s 
quark depending on whether the decay is a AS = or AS = —1 process. The operators Of are defined as 

0\ f = q a ^Lfpf^Lb a , O q 2f = q lfl Lff^Lb, 

Oy = q 7lM Lbq\L(R)q' , 0| >6 = q a ^Lbp<?^L{R)q' a , (2) 

3 3 
°7,9 = ^Ql^Lbeq'q'^Ri^q' , O| 10 = -q a -f^Lb p e q i^^R{L)q' a , 

where R(L) = 1 ± 75, and q' is summed over all flavors except t. Oi/,2/ ar e the tree level and QCD 
corrected operators. O^-q are the strong gluon induced penguin operators, and operators O7_io are due 
to 7 and Z exchange (electroweak penguins), and "box" diagrams at loop level. The Wilson coefficients c{ 
are defined at the scale [i ~ m& and have been evaluated to next-to-leading order in QCD. The c* are the 
regularization scheme independent values obtained in Ref. 0. We give the non-zero c{ below for nit = 176 
GeV, a s (mz) = 0.117, and fi = nib = 5 GeV, 

c x = -0.307 , c 2 = 1.147 , 4 = 0-017 , c\ = -0.037 , c{ = 0.010 , 4 = -0.045 , 
c\ = -1.24 x 10" 5 , 4 = 3.77 x 10" 4 , 4 = -0.010 , c{ = 2.06 x 10~ 3 , 
4$ = -c u A f 6 /N c = P^/N c , 4; 9 c = Pr, <*fa = (3) 

where N c is the number of color. The leading contributions to Pg e are given by: Pg = (|^)c2( 1 ^ + 
G(nii,^,q 2 )) and P % e = (^ L )(N c ci + c 2 )( 1 § + G(nii,^,q 2 )). The function G(m,/x,g 2 ) is given by 

G(m,iJ,,q 2 ) = 4/ x(l - a;) In ^ —dec. (4) 

jo M 

All the above coefficients are obtained up to one loop order in electroweak interactions. The momentum 
q is the momentum carried by the virtual gluon in the penguin diagram. When q 2 > 4m 2 , G(m, /j,,q 2 ) 
becomes imaginary. In our calculation, we use m u = 5 MeV, = 7 MeV, m s = 200 MeV, m c = 1.35 
GeV [||, [|. For T — » BX S , the operators Off and Orj/ do not contribute and the gluon momentum is fixed 
at q 2 = M\. 

A similar expression for the standard model contribution to the flavor changing decays J/ip can be 
written down. 
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To the standard model contribution we add higher dimensional four quark operators generated by 
physics beyond the standard model. In this paper, we consider the four quark operators with two scalar 
currents. 

Ken, = ^(l- 7 5 )^(l+7 5 )^+||s(l+7 5 )W(l-7 5 )b + /i.c. (5) 

The four quark operators in L new are the product of two scalar currents. In Eq. (5) A represents the 
new physics scale and R\ and R2 are two free parameters which describe the strength of the contribution 
of the underlying new physics to the effective operators. In our analysis we will only keep dimension six 
operators suppressed by 1/A 2 and neglect all higher dimension operators. 

Using a Fierz transformation one can express the scalar-scalar combination in terms of vector-vector 
combination. For instance we can write 

^(l-7 5 )^( 1 + 7 5 )^ = -^ a 7 M (l + 7 5 )6^(l- 7 5 )6 Q 

1 s Q 7 M (l+7 5 )6 a 6^(l-7 5 )6 /3 



2N C 



- ^(l-^J&^Wll-T 5 )^ (6) 

where T a are the SU(3) color matrices with the normalization Tr[T a T b ] = 5 a b/2 and N c is the number of 
colors. In the quarkonium system the leading component in the Fock space expansion involves the quark- 
antiquark pair being in a 3 5i state, probed by the operator S7 M (1 + 7 5 )Wry /i (l — 7 5 )6. Note the general 



Fock space expansion of quarkonium in NRQCD is |1C] 



|^q> = OWQQfS^ + OWlQQfpPlg) 

+ 0(v 2 )\QQ [ 3 S^ 8) ] gg) + 0(v 2 )\QQ [ l S^] g) + 0(v 2 )\QQ [ 3 D^ 8) ] gg) + • • • , (7) 

where v is the velocity of the constituents in the quarkonium and g represents a dynamical gluon, i.e. 
one whose effects cannot be incorporated into an instantaneous potential and whose typical momentum 
is tuqv 2 . The low energy hadronization of the leading component in the Fock space expansion of the 
quarkonium takes place at 0(v 3 ). As to the other Fock states notice that the | QQ [ 3 -Pj 8 ^] <?), configuration 



arises when the predominant state radiates a soft dynamical gluon. Such a process is mediated principally 
by the electric dipole operator, for which the selection rule is L' = L ± 1, S' = S, and which involves 
a single power of heavy quark three-momentum. Thus, the coefficient associated with this state is of 
order v. The electric dipole emission of yet another gluon involves a change from the P-wave state to the 
S- and D-wave states | QQ [ 3 Si ]gg), \ QQ), [ 3 Dj 8 ']gg] and so the associated coefficients are of order 
v 2 . Finally, the coefficient of the state | QQ \^Sq ] g) results from fluctuations into this spin-singlet state 
from the predominant spin-triplet state with the emission of a soft gluon via a spin-flipping magnetic 
dipole transition. Such transitions involve the gluon three-momentum (~ m,QV 2 ) rather than the heavy 
quark three-momentum (~ rriQv), and therefore the associated coefficient is of order v 2 . The low energy 
hadronization of these component in the Fock space expansion of the quarkonium takes place at 0(v 7 ). 
Also for the P wave, an additional factor of v comes from the derivative of the wavefunction. 

Before concluding this sections, we point out that besides operator in Eq. (5) there are additional four 
quark operators in the effective lagrangian|ll], such as those with vector- vector current structure, which 



contribute also to the processes we consider. We focus on operators in Eq. (5) because, as mentioned 
earlier, they can be generated by the exchanges of the new scalar bosons in models we consider below in 
section 5. 



3 Matrix Elements for T — > BX 



•S 



We proceed to calculate the matrix elements of the form < BX s \H e fj\T > which represents the process 
T — > BX S and where H e ff has been described above. The effective Hamiltonian consists of operators 
with a current x current structure. Pairs of such operators can be expressed in terms of color singlet and 
color octet structures. The factorization formalism based on NRQCD M, which allows a systematic and 
consistent probe of the complete quarkonium Fock space, can then be used to calculate the T decay rate. 

The matrix element of T —* BX S decay, can be expressed as, 
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where 



with 



M = < BX s \s^(l -7 5 )6|0 >< 0\b^b\T > 
v2 

+ < BX S \ «y (1 + 7 5 ) b I >< 0|6 7m 6|T > 
v2 

+ < «y (1 - 7 5 )T a 6 1 >< 0|6 7M T a 6|T > 
V 2 

+ < «y (1 + 7 5 )T a 6 | >< 0|6 7M T a 6|T > 
v2 

+ ^U 8 < BX S \ «y (1 - 7 5 )T a 6 | >< 0|6 7 m7 5 ^^|T > 
v2 

+ < 5X S | ay (1 + 7 5 )T a 6 | X 0|6 7At7 5 r a 6|T > 

+ ^=V 8 < BX s \s(l + 7 5 )T a 6 1 >< 0|6(1 - 7 5 )T a 6|T > 
v2 

+ ^y 8 ' < BX s \s(l - 7 5 )T a 6|0 >< 0|fr(l + 7 5 )T a 6|T > 
v2 



(8) 





= Wutd + Wi new 




= w{ std + w[ new 




= Wsstd + Wsnew 


wi 


= wi td + w^ new 


U s 


= Usstd + Usnew 


U's 


= u' 8std + u' Snew 


Vs 


= Vsstd + Vsnew 


V s ' 


= Vlstd + ^8new 



(9) 



{A 3 + A 4 - i(A 9 + A 10 )}(1 + + (As 



A fi 



iV 2 



(A 7 + 



= 



Wsstd 



[2(A 3 + A 4 + A 6 ) - (A 8 + A 9 + Aio)] 
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wL td = o 



U 8s td = [2(-i4g - A 4 + A 6 ) - (A 8 - A 9 



and 



We have defined 



with 



ULtd = o 



Vsstd = [-4A 5 + 2A 7 ] 

vLd = o (io) 



Wi new 



W 8 , 



1 V2R 2 
'2N C G F A 2 

1 v/2-Ri 
'2N C G F J? 
V2R2 
G F A 2 

A 2 



C/8new ~ G F A 2 



8new 



V2Rl 

G F A 2 



^8 = (11) 



^ = - E ( 12 ) 

(jr=M, c,t 



V q = V* s V qb (13) 



Similar expressions can be written for the matrix elements describing the J/tp decay. 
To calculate the decay rate we use the parton model to write the process T — > BX S as T(P) 
6(^1)5(^2) The squared matrix element is then given by 

C3; 



I MP = 2MZi 



< T|Oi( 3 Si)|T > (|H^i| 2 + \W[ 2 \) + (< T|0 8 ( 3 ^) + 2 ° s( 2 Pl) |T >)(|^ 8 | 2 + |wf 



m 2 h 
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+ 6MZ 2 [(< TIOg^lT > (|^s| 2 + K 2 \) 
+ 6MZ 3 [(< T\0 8 CS Q )\T > (\V 8 \ 2 + |vf |)" 



(14) 



where 



Z 2 



Zi 



= 8 -pi ■ p 2 + 



= 8 pi ■ p 2 + 



2pi • Pp 2 • P~ 
M 2 
2 Pl -Pp 2 -P 
f M 2 



Z 3 



= 8 [pi • P2] 



(15) 



with M being the quarkonium mass. The matrix elements of the various color singlet and color octet 



Along with the CP violating phases present in the standard model contribution there can be additional 
phases from the new contribution. We can then construct the CP violating rate asymmetry as 



4 Low Energy Constraints 

The lagrangian L new generates, at one loop level, the effective S67*, sbg* ,sbZ vertices as shown in Fig. 1, 
where 7* and g* indicate an off shell photon and a gluon.. Similar vertices involving c — > u transitions are 
generated in the charmonium sector also. These vertices, with a 7 and Z, will contribute to b — > sl + l~ 
and c — > ul + l~ . Note there is no contribution to b — ► 57. The vertex b — ► s<7* can give rise to the process 
6 — > sgg which will contribute to non-leptonic B decays. We expect the constraints from b — > to be 

better than from non-leptonic -B decays because of the theoretical uncertainties in calculating non-leptonic 
decays. The effective sbj*, sbg*,sbZ lagrangian can be written as 




a C p = 



r(T BX S ) - T(T BX S ) 
r(T ^ + T(T -> BX S ) 



(16) 



5L 3&7 * = -^[R+Cl + R.C^b 



(17) 
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Y*, g*, Z 




where 



R + 
R- 

C+ 



Fi gure 1: Effective sb^f*, sbg* 7 sbZ vertices generated by L Tl 



R\ + i?2 

2 

R\ — i?2 



2 



^2 ^ (^2 ) [ 8 9 M 7 • 97 5 ^ - 1) + 8 7 ^ 7 5 g 2 x(l - x) 



(18) 



with 



= ml - q 2 x(l - x) 



where ej, is the b quark electric charge and q is the photon momentum. The effective sbg* lagrangian can be 
obtained by replacing eb by g s , the strong coupling constant. The effective sbZ lagrangian can be written 
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as 



where 



Fi- 
F 2 - 



with 



Sir. 



bZ 



9 - 



2c w A 2 



s [R+(g L F 1+ + g R F 2+ ) + R-(g L F^ + g R F 2 -)] 



8<f 7 • qx(x - 1)(1 + 7 5 ) + 8^q 2 x(l - x)(l + j 5 ) - tof'fm* 



1 



167T 2 Jo 

1 rt 



dxlog —5 8g^7 • gx(:r - 1)(1 + 7 5 )87^g 2 a;(l - x)(l + 7 5 ) - &fm\ 



A 2 ' 

jf 1 dxlog (^) [8^7 • - x)(l - ^)8^q 2 x(x - 1)(1 - 7 5 ) + 8 7 "m 2 



f 7 • gx(:r - 1)(1 - 7 5 )87^g 2 x(l - x)(l - 7 5 ) + 87^7 



5 m 2 



9l 
9r 



1 1 2 

"2 + 3 S ™ 



(19) 



(20) 



For b — > sl + l~ the terms in the equations above do not contribute if we neglect the lepton masses. 
Furthermore the contribution from the Z exchange is suppressed with respect to the 7 exchange by factor 
of q 2 : /M§ and so we do not include the Z contribution. The additional contribution to the effective 
Hamiltonian for b — * sl + l~~ can be written as 

A 2 " 



SH, 



b->sl+l- 



^L-^fidxteil-z) log . - 2 



R^bLljJ + R2vfb R liJ 



21 



which has to be added to the standard model contribution ||. Similar results can also be written for the 
charm sector. 



5 Models 



In this section we look at various models that can give rise to I new given in Eq. 5. As a first example we 
consider a recent version of top color models [13]. In such models the top quark participates in a new strong 
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interaction which is broken at some high energy scale A. The strong interaction, though not confining, 
leads to the formation of a top condensate < tit R > resulting in a large dynamical mass for the top quark. 
The scale A is chosen to be of the order of a TeV to avoid naturalness problem which implies that the 
electroweak symmetry cannot be broken solely by the top condensate. In the low energy sector of the 
theory, scalar bound states are formed that couple strongly to the b quark |13|, PJ 

L b = -P-b L {H + iA°)b R + h.c (22) 

JnV2, 

where ~ 50 GeV is the top pion decay constant. On integrating out the Higgs fields H and A we have 
an effective four fermion operator 

2 

Yfi 

L eff = -72Z^b L b R b R b L (23) 
Since the b quark in (22) is in the weak-eigenstate, L e ff in (22) will induce flavor changing neutral current 



(FCNC) four quark operators in Eq. (5) after diagonalizing the quark mass matrix [14], with coefficients, 

Rl = A f2 *2 \ D Lbb\ 2 D Rbb D* Rbs 

R * = \ 1 P^r\D R ob\ 2 D Lhb Di hs (24) 

where Dl and D R are the mixing matrices in the left and the right handed down sector. In the charm 
sector similar interactions can arise due to the strong couplings of the top quark to top pions. The effective 
operators generated by integrating out the top-pions are similar to Eq. (5) with replacement of b by c and 



s by u. In topcolor II models [14, 15], where there can be strong top-pion couplings of the top with the 
charm quark, we have 

1 m 2 



Rl ~ Jj2^2\ U Lcc\ 2 U Rt cU Rtu 

^ = \^\Um c \ 2 U Lcu U* Lcc (25) 

J TV 7T 



In supersymmetric standard models without R parity, the most general superpotential of the MSSM, 
consistent with SU(3) x SU(2) x U(l) gauge symmetry and supersymmetry, can be written as 

W = W R + W r , (26) 
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where Wr is the -R-parity conserving part while Wj^ violates the -R-parity. They are given by 

W R = hijLiHiEf + ^QiHiDl + hZjQiHxUf, (27) 
W# = KjkLiLjEl + X' ijk LiQjDf, + X ijk UiDjD k + ^Lji^. (28) 

Here Lj(Qj) and Ei(Ui,Di) are the left-handed lepton (quark) doublet and lepton (quark) singlet chiral 
superfields, with i,j,k being generation indices and c denoting a charge conjugate field. Hi 2 are the 
chiral superfields representing the two Higgs doublets. In the -R-parity violating superpotential above, 
the A and A' couplings violate lepton-number conservation, while the A" couplings violate baryon-number 
conservation. Xijk is antisymmetric in the first two indices and X"j k is antisymmetric in the last two indices. 
While it is theoretically possible to have both baryon-number and lepton-number violating terms in the 
lagrangian, the non-observation of proton decay imposes very stringent conditions on their simultaneous 
presence We, therefore, assume the existence of either L- violating couplings or l?-violating couplings, 
but not the coexistence of both. We calculate the effects of both types of couplings. 

In terms of the four-component Dirac notation, the lagrangian involving the A' and A" couplings is 
given by 



Cy — -X' ijk 



444 + 444 + (4)*K) C 4 
-444 - 444 - (4)*(4) c 4] + ( 2 Q) 
c x „ = -a^[4(4) c 4 + 4(4) c < + 4(4) c 4]+^-c. (30) 

The terms proportional to A are not relevant to our present discussion and will not be considered here. 
The exchange of sneutrinos with the A' coupling will generate L new for T — > BX S with 

r> x v A a2 A m 

K\ — -Zji o — 

4 mf 7 

R 2 = (31) 
For the case of J ftp — > DX U the operators in L new cannot be generated at tree level. 
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Another model of interest is an extension of the SM with additional scalar SU(2) doublets, the simplest 
of these would be the two Higgs doublet model (2HDM). In general, when the quarks couple to more 
than one scalar doublet, there are inevitably FCNC couplings to the neutral scalars. When the up- 
type quarks and the down-type quarks are allowed simultaneously to couple to more than one scalar 
doublet, the diagonalization of the up-type and down-type mass matrices does not automatically ensure the 
diagonalization of the couplings with each single scalar doublet. Frequently, as in the Weinberg model for 



CP violation [17] or in Supersymmetry, the 2HDM scalar potential and Yukawa lagrangian are constrained 
by a ad hoc discrete symmetry fig}] , whose only role is to protect the model from FCNC's at the tree level. 
Let us consider a Yukawa lagrangian of the form 

£y } = riYjQi^iU^R + r]ijQi,L<j>iDj i R + ^Qi^U jtR + ^Qi,L^D j:R + h.c. (32) 

where <f>i, for i = 1,2, are the two scalar doublets of a 2HDM, while vfj D an d ^fj D are the non-diagonal 
matrices of the Yukawa couplings. 

When no discrete symmetry is imposed then both up-type and down-type quarks can have FC couplings 



19]. Such models were called Class A in (2(J to be contrasted with models with a forced absence of FCNC, 
called Class B. 

In the notation and basis of Ref[^l|] the flavor changing part of the lagrangian can be written as 

£y"c = £YjQi,LfaU jtR + i^Qi^hD^R + h.c. (33) 

where Qi,L, Uj : R, and Dj : r denote now the quark mass eigenstates and are the rotated couplings, 
in general not diagonal. If we define V^'^ to be the rotation matrices acting on the up- and down-type 
quarks, with left or right chirality respectively, then the neutral FC couplings will be 

Wutral \ L ) •£ ' V R ■ I 34 ) 
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On the other hand for the charged FC couplings we will have 

^charged ^neutral ' ^CKM 

^charged = ^CKM ■ luteal (35) 

where Vckm denotes the Cabibbo-Kobayashi-Maskawa matrix. 

The phenomenology for the 2HDM, for the quarkonium processes under study, is not very different 
from the other two models considered above and so we will concentrate mainly on the top color models 
and supersymmetry models with R-parity violation. 

6 Results and Discussion of Theoretical Uncertainties 



Table 1: Input parameters used in our calculations. 



NRQCD matrix elements 


Value 


(Of( 3 5 1 ))«3(^|Oi( 3 5 1 )|^) 
(TIO^SOIT) 

(otesi)) 

(T|O g ( 3 Si)|T) 

(oieso)) « (oje Po ))/mi 


0.73 GeV 3 
2.3 GeV 3 
0.014 GeV 3 
10- 2 GeV 3 
5 x 10~ 4 GeV 3 
7 x 10~ 3 GeV 3 



In this section we discuss the results of our calculations. First let us look at the T decays. The inputs to 
our calculation are the various well known NRQCD matrix elements given in table. 1 |f22| . The standard 
model contribution to the branching ratio is 5.2 x 10~ n from the penguin induced b — ► s transition. The 
process T — * BX S can also have a contribution in the standard model from tree level processes. The 
effective Hamiltonian , suppressing the Dirac structure of the currents, 

G F _ _ 

H w = —=V ub V* s [ai(ub)(su) + a2(sb)(uu)] (36) 
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where a\ and are the QCD coefficients can generate the process T — > B + K . We can estimate the 
branching ratio for this process as 

BR[T -> fl+K"] « |^| 2 Si2[T -> 
Vc6 

Using BR[T -> B+2f~] calculated in Ref@ one obtains Bi?[T -> fl+K"-] ~ 1.5 x 10" 14 . For a rough 
estimate of BR[T -> £+X s ] we can scale RR[T -> B+if - ] by the factor D°X}/BR[B -> Dtt]. 

The measured value of BR[B — > D°X] || includes Z?° coming from the decay of -D * and From the 

spin phase factors AR[B -> L>*X] ~ WR[B -» DX\. Hence BR[B -» D°X}/BR[B -> Dtt] ~ 20 leading 
to i?i?[Y -> 5+X s ] ~ 3 x 10~ 13 . 

So far we have not considered Ri and i?2 in L new . In our model independent analysis we vary Ri/A 2 , 
R2/X 2 one at a time and the use the constraint from measurements of b — > se + e~ and 6 — ► sfi + n~ . We 
identify A with the masses of the exchange particles which we take to be between 100 — 200 GeV. We also 
take the cut-off for the integral in Eqs. 17-20 as 200 GeV. The allowed values of Ri/A 2 , R2/A 2 are then 
used to calculate T — ► BX S The constraint from b — > sl + l~ gives 

|i?i, 2 |/A 2 < (6 - 9) x lQ-*(l/GeV) 2 

Using the upper bounds on |i?i2|/A 2 we find the branching ratio for the process T(15) — > BX S to be 
between (1 — 2) x 10~ 6 . Branching ratios of similar order are also obtained for T(25) and T(3S). For 
T(45) the branching ratio is smaller by a factor of 100 because of the larger width of Y(4f>) which decays 
predominantly to two B mesons. 

Turning now to models, we find for the top color model from Eq. (24) we can write 

DRbs ~ ^ml\D Lbb \*D Rbb (37) 

We can identify A = mu and use the constraint from b — > se + e~ for a typical value of |i?i|/A 2 ~ 
6 x lO^l/GeU) 2 . Assuming \D Lbb \ ra \D Rbb \ 1, and / ff = 50GeV we obtain \D Rbs \ ~ 2m 2 H x 10~ 6 . 
With typical values of mjj ~ 100 — 200 GeV we get l-D^bsl ~ 0.02 — 0.08. Similar values have been 
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obtained for | Df>b s | in Ref |14| by considering the contributions of the charged higgs and top-pion to b — > 
S7. A similar exercise can be carried out with |Z?ife s |. Note that B s mixing probes the combination 
D^gD RbbD* Rbs D ibb and so by either choosing R\ ~ or R2 ~ we can satisfy the constraint on i? s mixing 
by choosing the appropriate mixing elements to be small. Note that in top color models we can have 
operators s(l — 7 5 )M(1 + 7 5 )d and s(l + 7 5 )6c?(l — 7 5 )d that can contribute to T — > Bsd — > BX S after Fierz 
reordering. However these operators will be suppressed by form factor effects and also from mixing effects. 
We have checked that the contribution to T — > BX S from these operators are much suppressed relative to 
the contribution of the operators in L new . We will therefore not consider the the above operators in our 
analysis. 

Turning to R-parity violating susy we first collect the constraints on the relevant couplings. The upper 
limits of the L-violating couplings for the squark mass of 100 GeV are given by 

L,2), (38) 

(39) 
(40) 
(41) 
(42) 

The first set of constraints in Eq. (38) come from the decay K — ► m>v with FCNC processes in the down 



\l 1 

^kij 1 


< 


0.012, (k,j = l,2,3;i 


A 'l3jl 


< 


0.16, (j = 1,2), 


A 'l33 


< 


0.001, 


A 23i ] 


< 


0.16, (j = 1,2,3), 


A 33i 1 


< 


0.26, (j = 1,2,3), 



quark sector [24]. The set of constraints in Eq. (39) and Eq. (41) are obtained from the semileptonic 
decays of l?-meson [p5f| . The constraint, on the coupling A' 133 in Eq. (40) is obtained from the Majorana 
mass that the coupling can generate for the electron type neutrino pq| . The last set of limits in Eq. (42) 
are derived from the leptonic decay modes of the Z [^7]] . Assuming all the couplings to be positive we find 
the branching ratio for T — > BX S to be around 2 x 10~ 6 for = lOOGeV. 

Turning next to J ftp — > DX U , we first make an estimate for this process in the standard model. Since 
the penguin c — > u transition is small in the standard model we neglect its contribution. As in the case for 
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the T system , for a rough estimate, can write 

BR[J/t/j D°X U ] ~ BR[J/ip D°-k°}BR[D K~X}/BR[D -> if _ 7r + ] 

We obtain BR[J/tp — ► D°7r°] from j|] and keeping in mind that — > contains contributions 

from states decaying to if - we obtain BR[J '/ip — > _D°X U ] ~ 10~ 10 . A similar exercise gives BR[J/ip — > 
D+X u ] ~ 10- 9 . 

Considering new physics effects we can constrain i?i and i?2 from c — > ul + l~ . We get an estimate of 
the constraint on c — ► iie + e~ by adding up the exclusive modes 

5i?[L> -» ue+e"] > BR[D -» (tt° + 7/ + p° + w)e + e~] 

From c — ► ul + l~ one obtains 

|^i, 2 |/A 2 < 3.7 x KT 4 (l/GeF) 2 

We find the branching fraction for the process J/t/j — > DA U using the constraint from c — ► ul + l~ can be 
(3 - 4) x 1(T 5 

In top color models taking R\ and i?2 one at a time, one obtains 

2.1 x 10 3 , 



* |2 



mi 



■\\Ulcc\ U Rtc U* Rtu 



or 



4 — 1 1 C^ffic 1 2 U L ccU* Lcu 1 2 

as the branching fraction for J ftp — > DX U . For mj between 100 — 200 GeV this rate can be between 
(0.1-2.0) x 10" 5 if all the mixing angles are ~ 1. It has been shown in Ref[^8| that our choices for and 
m% gives unacceptably large corrections to Z — > bb from one loop contribution of the top pions. However 
in a strongly coupled theory higher loop terms can have significant contributions. Nonetheless if we change 
fjr to ~ 100 GeV for better agreement with Z — » bb data then the effect in J/tp — * DX U is reduced by 
a factor of 16. As in the case of the T system we can satisfy the constraint from D mixing by choosing 
Ri ~ or R 2 ~ 0. 
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For R parity violating susy, contribution to J/ijj — > DX U can only occur at loop level, with both the 
A' or A" contributing, through the box diagram and so is suppressed. However in the general 2HDM, from 
Eq. 33 , the operator uccc can be generated by the tree level exchange of the field 02- The contribution 
to J/tp — > DX U will be proportional to the combination of couplings R2 = Ci2^22* an d -Ri = ?2i*^22- Note 
the D° — Z>° mixing probes ^12^21*- So we can satisfy the constraint on D° — D° mixing by choosing either 
Ri ~ or R2 ~ 0. One then obtains 

14 \£U £U*\2 
"IS12S2 



m 



4~|S12S22 
H 



or 



14 \£U*tU\2 
~ls21 S22l 
H 



rn 



as the branching ratio for J/ip — ► DX U . For m# between 100 — 200 GeV this rate can be between 
(0.1 — 1.4) x 10~ 7 if all the couplings are ~ 1. In a strongly interacting theory these couplings can be > 1 
as in the example of top color models discussed above. 

We note in passing that the four quark operators we have considered can also give rise to mixing 
operators that are 1/A 4 suppressed. Taking one operator at a time one can generate the following operators 
that contribute to mixing 



We then have for B s mixing, 



A5 S = ^f Bs 2 r)BM Bs 5, (44) 



where 77 is the QCD correction factor and B and S are defined through 

< B° s \s(l ~ 7 5 M1 " 7 5 )l^° >= -\f Bs 2 BM Ba \ 



and 



3R 2 ml, fA 2 ^ 



5 2--'A-' A2 l0g \ n ,i 
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where i = 1, 2. A similar result can be written for D mixing. 

For the Upsilon system if we include the constraint from B s mixing generated by only the operators 
Oi : 2 then we obtain {Ri^l/A 2 > 1 — 2 x W~®(l/GeV) 2 . This is of the same order as the constraint obtained 
from b — * sl + l~ . In the case of J/tp, if we include the constraint from the D mixing generated by only 
the operators Oxp, then we obtain |i?i 2|/A 2 ~ 10~ 6 (l/GeV) 2 . This will lower the branching fraction for 
J/i/j — ► DX U to ~ 10~ 9 . However,to evaluate consistently the effects of these operators at order of (1/A 2 ) 2 
would require the addition of other operators with dimension < 8. As noted earlier we restrict our analysis 
to only dimension six operators in the effective lagrangian, thus the conservative result for J/i/j — > DX U is 
of order of 10 -5 , which, as we shown above, lies also in the region predicted by the TopColor and 2HDM. 

Direct CP violation is possible in these decays through the interference of the standard model and new 
physics contribution. The CP conserving phase is generated at the quark level from the penguin diagrams. 
However the standard model contribution is small and so the CP asymmetry acp is also small with a 
typical value of 0.1% for T — ► BX S . 

In summary we have calculated branching ratios for the flavor changing processes T — ► BX S and 



J/ij) — > DX U . In a model independent description of new physics [29], constrained by low energy data from 
b — * sl + l~ and c — * ul + l~ , we found branching fractions for these processes can be ~ 10" 6 and ~ 10~ 5 
for T and J/ijj decays respectively. We also discussed several models of new physics that can allow these 
processes to occur with branching ratios that maybe measurable in the next round of experiments. 
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